Improving entanglement concentration of Gaussian states by local displacements 
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We investigate entanglement concentration of continuous- variable Gaussian states by local single- 
photon subtractions combined with local Gaussian operations. We first analyze the local squeezing- 
enhanced entanglement concentration protocol proposed very recently by Zhang and van Loock [e- 
print: arXiv:1103.4500 (2011)] and discuss the mechanism by which local squeezing before photon 
subtraction helps to increase the entanglement of the output state of the protocol. We next show 
that a similar entanglement improvement can be achieved by using local coherent displacements 
instead of single-mode squeezing. 

PACS numbers: 03.67.Mn,03.67.Hk, 42.50.Dv 
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I. INTRODUCTION 



Distillation of continuous- variable entanglement is fun- 
damentally limited by a no-go theorem which states that 
experimentally accessible entangled Gaussian states can- 
not be distilled by local Gaussian operations and clas- 
sical communication [lH3|. I n order to increase entan- 
glement of a Gaussian state such as two-mode squeezed 
vacuum by local operations, one thus has to resort to 
non-Gaussian operations [4J]. A prominent example of 
such an operation is the photon subtraction, that has 
been proposed for continuous- variable entanglement con- 
centration in a seminal paper by Opatrny et al. [5|- 
Since then, the photon subtraction has been recognized 
as a crucial tool for implementing many tasks in opti- 
cal continuous- variable quantum information processing. 
Besides entanglement concentration |6h10||. this includes 
generation of Schrodinger cat-like states by subtracting 
photons from single- mode squeezed vacuum state flll - tl4| , 
generation of arbitrary single-mode states of light beams 
[15l 1 161 ] , implementation of quantum gates for qubits en- 
coded into superpositions of two coherent states [17j , and 
probabilistic noiseless amplification of coherent states of 
light HH23- 

Although Gaussian operations alone are useless for dis- 
tillation of Gaussian entanglement, it has been recently 
shown by Zhang and van Loock [25[ that the perfor- 
mance of entanglement concentration via local single- 
photon subtraction can be enhanced by suitable local 
squeezing operations. In this paper we further inves- 
tigate this intriguing concept and clarify the origin of 
enhancement of entanglement concentration due to the 
local pre-squeezing. We next show that the same ef- 
fect can be achieved also by local displacements, which 
are much easier and experimentally more feasible opera- 
tions than squeezing. Indeed, the combination of photon 
subtraction and coherent displacement has been recently 
successfully experimentally utilized to generate arbitrary 
superpositions of squeezed vacuum and squeezed single- 
photon states from squeezed vacuum J16fl. Coherent dis- 
placement can be implemented by mixing the signal beam 
with a weak coherent laser beam on a highly unbalanced 
beam splitter. In contrast, squeezing is much more exper- 



imentally challenging and resource demanding operation. 
Squeezing in cavity-based optical parametric amplifier 
would require efficient injection of the pulsed signal beam 
in the cavity as well as efficient retrieval of the squeezed 
beam which is a nontrivial task. Squeezing based on a 
single passage through a nonlinear crystal pumped by a 
strong pulsed laser still requires careful modematching 
and may suffer from noise due to spontaneous paramet- 
ric down-conversion producing pairs of photons in many 
modes. Alternatively, one can employ an off-line gen- 
erated auxiliary squeezed vacuum state combined with 
interference on an unbalanced beam splitter, homodyne 
detection and feedforward [26|, [27| ■ This squeezing tech- 
nique was successfully experimentally tested, however, 
the operation is imperfect and inevitably adds some ex- 
tra noise due to finite squeezing of the auxiliary state. 

The rest of the paper is structured as follows. In Sec- 
tion II we present a simple analysis of the local squeezing- 
enhanced entanglement concentration which reveals the 
origin of the entanglement enhancement in this scheme. 
In Section III we describe the displacement-enhnaced en- 
tanglement concentration procedure and we analytically 
determine the optimal displacement and success proba- 
bility scaling in the limit of pure weakly squeezed input 
states, where the scheme achieves the maximum entan- 
glement enhancement. We also numerically calculate the 
optimal local displacements maximizing entanglement of 
the distilled state for both pure and mixed Gaussian two- 
mode squeezed states. Finally, Section IV contains brief 
summary and conclusions. 



II. SQUEEZING-ENHANCED 
ENTANGLEMENT CONCENTRATION 

The considered entanglement concentration scheme is 
depicted in Fig. 1(a). The entangled two- mode squeezed 
state is distributed from a source via some quantum chan- 
nel to the two recipients traditionally called Alice and 
Bob who attempt to conditionally increase the entangle- 
ment of the shared state by local photon subtraction, 
which is a specific example of a quantum filter. The 
whole scheme is thus a typical example of a Procrustean 
entanglement concentration protocol [28]. In practice, 
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pre-squeezing can be useful, we apply an inverse unitary 
squeezing operation S\(s) (g) S' B (s) to the concentrated 
state which does not change its entanglement, see Fig. 
1(c). The overall filtering operation Fab then reads 
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FIG. 1: (a) Entanglement concentration by single-photon sub- 
traction. BSa and BSb are highly unbalanced beam split- 
ters, pi n is the shared (generally mixed) entangled state, and 
Da and Ds denote single-photon detectors, (b) Entangle- 
ment concentration improved by local single-mode suqeezing 
S. The state can be un-squeezed after concentration by in- 
verse local operations §' without changing its entanglement 
(c). 



the photon subtraction is accomplished by means of a 
highly unbalanced beam splitter that reflects a tiny part 
of the incoming beam on a single-photon detector APD. 
A successful photon subtraction is heralded by a click 
of the detector. Alice and Bob exchange their measure- 
ment results via classical communication line and the en- 
tanglement concentration succeeds if they both register 
detector clicks. In the limit of very low beam splitter 
reflectance, the resulting operation can be described by 
action of annihilation operators a and b on modes A and 
B, respectively, even if the detectors have non-unit detec- 
tion efficiency and cannot distinguish number of photons. 
It is instructive to first investigate concentration of a 
pure two-mode squeezed vacuum state 



|V^> = V / l-A 2 ^A"|nn), 



(1) 



71=0 



where \n) denotes Fock state, \nn) is a short-hand nota- 
tion for |n)|n), and A = tanhr, where r stands for the 
two-mode squeezing constant. The non-normalized state 
after local photon subtractions can be expressed as 



Fab = S\{s)aS A {s) <g> Sl(s)bS B (i 



(3) 



The squeezing transforms the annihilation operators into 
linear combination of annihilation and creation opera- 
tors, and we have 

Fab = [acosh(s) + a^ sinh(s)] <g> [6cosh(s) + 6' sinh(s)]. 

(4) 
We can thus see that the local pre-squeezing effectively 
converts the photon subtraction into a coherent combi- 
nation of photon subtraction and photon addition. 

Let us now investigate what happens in the limit of 
weak two-mode squeezing, A< 1. We can write 



FabW) = Fab (|00) + A|ll» + 0(A 2 ), 



(5) 



where 0(A 2 ) indicates small terms of the order of A 2 . 
After some algebra we get 

F AB \ip) = Acosh 2 (s)|00)+smh 2 (s)|ll)+2Asinh 2 (s)|22) 
A 



V2 



sinh(2s)(|02) + |20)) + O(A 2 



If we set tanh(s) = VA, then we obtain 
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Note that the amplitude of the states |02) and |20) is 
suppressed by a factor of y/\ with respect to the am- 
plitude of the dominant term |00) + |11). In the limit 
of very small A the result of entanglement concentration 
is a maximally entangled state of two qubits hence the 
protocol can extract one e-bit of entanglement from ar- 
bitrarily weakly squeezed initial state. The success prob- 
ability of the protocol scales as P S ucc oc A 2 for A <C 1. 
The squeezing enhanced protocol should be contrasted 
with the ordinary concentration by local photon subtrac- 
tion that in the limit A < 1 provides only very weakly 
entangled state |00) + 2A|11). 



III. DISPLACEMENT-ENHANCED 

ENTANGLEMENT CONCENTRATION 
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Alice and Bob may amend the entanglement concentra- 
tion protocol by local squeezing of modes A nad B prior 
to the photon subtraction [25], c.f. Fig. 1(b). The fil- 
tering operation then reads a (3 bSA{s) <S> Sb(s), where 
§m is a single-mode squeezing operation on mode M 
with squeezing constant s. In order to clarify why the 



We next show that similar improvement in the amount 
of extracted entanglement can be obtained if we replace 
the local squeezing operations with local displacements 
Da{o) and Db(P), see Fig. 2(a). Similarly as above, 
we assume that the displacement is undone after photon 
subtraction by inverse displacements D A (a) and D' B (f3), 
c.f. Fig. 2(b). The resulting local quantum filter reads 



Gab = D A (a)aD A {a) <g> D B (fi)bb B {P). 



(7) 
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FIG. 2: Entanglement concentration scheme combining lo- 
cal displacements and photon subtraction (a). Applying lo- 
cal unitary inverse displacement operations on the state after 
concentration does not change its entanglement (b). Instead 
of displacing the signal modes one can equivalently inject 
weak coherent states into the auxiliary ports of the unbal- 
anced beam splitters (c) or displace the beams impinging on 
single- photon detectors (d). 



The local displacements transform the photon subtrac- 
tion into a coherent superposition of photon subtraction 
and identity operation and we have 

G AB = (a + a)<E){b + l3). (8) 

Let us again consider the limit A <C 1. We obtain 

G AB \ip) = (a/3+A)|00)+A(a|10) +/3|01))+Aa/3|ll)+O(A 2 ) 

The apparently dominant vacuum term corresponds to a 
separable state, but it can be completely eliminated by 
a destructive interference provided that a/3 = — A. The 
maximal entanglement is then obtained if the displace- 
ments are real and symmetric, 

a = -0 = VA. (10) 

If we insert a and (3 given by Eq. (jTU)) into © we get 

GabW = A 3 / 2 (|10) - |01)) + 0(A 2 ). (11) 

We can see that in the limit Ac lwe again extract one e- 
bit of entanglement, although now in a slightly different 
form of a single photon split among the two modes A 
and B. The success probability of the scheme scales as 



fsucc oc A 3 for small A. In comparison with the pre- 
squeezing strategy the success rate is thus smaller by a 
factor of A. Nonetheless, the simplicity of displacement 
operation compensates for this reduced efficiency. 

Note that the local photon subtractions combined with 
displacements can be seen as a local version of the non- 
local coherent photon subtraction demonstrated by Our- 
joumtsev et al. 8]. There the two beams reflected off the 
beam splitters BS^ and BSb were combined on another 
balanced beam splitter and a click of an APD placed on 
one of its output ports heralded coherent photon sub- 
traction associated with the filter Gab,h\ — o, + b. Sur- 
prisingly, the present local scheme becomes equivalent 
to that non-local strategy in the limit of small A and 
for carefully tuned local displacements (TTU1) . Note also 
that the displacements can be equivalently performed on 
the ancilla input ports, e.g. by injecting weak coherent 
states there, see Fig. 2(c). Yet another possibility is to 
displace the reflected beam right before the single pho- 
ton detection, c.f. Fig. 2(d). This latter scheme has 
been recently successfully experimentally utilized to gen- 
erate qubit states formed by squeezed superpositions of 
vacuum and single-photon states [lq |. 

Let us now analyze the performance of the 
displacement-enhanced entanglement concentration 
scheme beyond the limit of A C 1. In our model 
we take into account a non-zero reflectance R of the 
beam-splitters BS^ and BSs and we assume on-off 
single-photon detectors described by a two-component 
POVM ri = |0)(0| and fti = I - fl , where I denotes 
the identity operator and II i corresponds to a click of the 
the detector. For the sake of simplicity we assume that 
the detectors have unit efficiency. Experiments usually 
involve highly unbalanced beam splitters with fi« 1. 
In this limit, the main effect of non-unit detection 
efficiency rj is the reduction of success rate by a factor 
of rj 2 . More specifically, the entanglement concentration 
scheme with inefficient detectors is equivalent to a 
scheme with perfect detectors and reduced reflectance 
of the beam splitters, R' = rjR/[l — (1 — rfjR], preceded 
by the transmission of the state through a lossy channel 
with transmittance T = 1 — (1 — rj)R. In the limit of 
flCl the latter can be neglected since TrjI. We allow 
for an arbitrary initial mixed two-mode Gaussian state 
Pi n as we shall investigate entanglement concentration 
of two-mode squeezed vacuum states transmitted over 
lossy channels. We denote by Ubs,ac and Ubs,bd the 
unitary operation describing mixing of two modes on 
the unbalanced beam splitters with reflectance R. The 
non-normalized density matrix of the concentrated state 
can be expressed as follows [25j [ 



Pout 



Tr, 



CD 



{licD ® Iab) U {pin,AB ® |00)(00| C d) U 1 

(12) 



where 

U = Ubs,ac <8> Ubs,bd Da{o) ® D B (/3) <E> I CD , (13) 
TLcd = III C ® III d, an d Trc£> denotes a partial trace 
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FIG. 3: Logarithmic negativity En of two-mode squeezed 
state after entanglement concentration. Contour plot in panel 
(a) shows dependence of En on the local displacements a and 
/3 and panel (b) represents the cut through panel (a) along 
the diagonal line /3 — — a. The parameters read A = 0.15 and 
# = 0.1. 



over modes C and D that are initially prepared in the 
vacuum state |00)c_d. 

We quantify the entanglement of the concentrated 
state pout by logarithmic negativity Erf, which is a com- 
putable measure of entanglement for non-Gaussian mixed 
states [29|. Note that /5 ou t is mixed even for pure input 
Pi n due to the inability of the APDs to resolve number of 
photons. Recall that 



E n (pab) = log 2 



\p T a a b\ 



(14) 



where Ta denotes partial transposition with respect to 
mode A and \\X\\i = TrvXtX. In order to calcu- 
late En we need to determine matrix elements of /5 ou t 
in Fock state basis. One option is to explore the fact 
that for a Gaussian input state p- m the output state /5 out 
can be expressed as a linear combination of four Gaus- 
sian states [lji [25|, [3(j. For a Gaussian state, the density 
matrix elements in Fock basis can be expressed in terms 
of multivariate Hermite polynomials because the Husimi 
Q-function of the state, which is a generating function 
of these matrix elements, has a Gaussian form. Here we 
follow a more straightforward route and directly evaluate 
Eq. ([T^jl in Fock basis truncated at Fock state n max = 10. 
For the considered values of A and other parameters, the 
errors due to the truncation are negligible. 

In Fig. 3 we plot the dependence of logarithmic nega- 
tivity of concentrated state on the local displacements 
a and /3 for initial pure two-mode squeezed vacuum 
with A = 0.15 and beam splitters with 10% reflectance, 
R = 0.1. We can see that the two peaks of maximum 
entanglement lie on the line a = — /3, as expected from 
the symmetry. Numerical calculations indicate that it 
suffices to consider real displacements, adding complex 
phase factors does not further increase the entanglement. 
In the rest of the paper we thus restrict ourselves to real 
a and always set j3 = — a. 

The dependence of the performance of entanglement 
concentration protocol on the two-mode squeezing con- 
stant A is illustrated in Fig. 4. The logarithmic negativ- 
ities of an initial pure two-mode squeezed vacuum state, 



a state after local single-photon subtractions, and a state 
after local single-photon subtractions combined with op- 
timal local displacements are shown in Fig. 4(a). We 
can see that the scheme with local displacements outper- 
forms the scheme involving only photon subtractions and 
it yields almost 1 e-bit of entanglement even in the limit 
A <C 1 . The slight reduction of En below 1 e-bit is caused 
by on-off single-photon detectors leading to mixed output 
state after entanglement concentration. The data plotted 
in Fig. 4(b) confirm that the success probability of the 
protocol scales as i? 2 A 3 (with displacements) and i? 2 A 2 
(without displacements) . The optimal displacement a op t 
follows the scaling a op t = VA for small A, c.f. Fig. 4(c), 
but for larger A it holds that a op t < \/A. 

Finally we investigate entanglement concentration of 
mixed states obtained by transmitting each mode of a 
pure two-mode squeezed vacuum state through a lossy 
channel with transmittance 1 — v. Results of numerical 
simulations plotted in Fig. 5 reveal that the local dis- 
placements help to increase the attainable entanglement 
even in the presence of losses although the achievable 
entanglement decreases with increasing channel losses v. 
Also the success probability of concentration and the op- 
timal displacement decrease with increasing losses, see 
Fig. 5(b) and 5(c). This effect can be attributed to sig- 
nal attenuation and reduction of mean number of photons 
due to transmission through a lossy channel. 



IV. CONCLUSIONS 

In summary we have investigated continuous- variable 
entanglement concentration schemes based on local pho- 
ton subtractions combined with local Gaussian opera- 
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FIG. 4: Performance of entanglement concentration in de- 
pendence on the two-mode squeezing constant A. In panel 
(a) we plot the logarithmic negativities En of three states: 
initial pure two-mode squeezed vacuum (solid line), a two- 
mode state after local photon subtractions (dotted line) and 
a two-mode state after local photon subtractions combined 
with optimal local displacements (dashed line). Panel (b) 
shows the success rates P SU cc of the concentration protocol 
with (dashed line) and without (dotted line) displacements. 
Panel (c) shows the optimal displacement a opt (solid black 
line) and the vA reference curve for comparison (gray dashed 
line). The results were obtained for R = 0.1. 
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FIG. 5: Entanglement concentration of a two-mode squeezed 
state transmitted over a lossy channel with loss factor v. 
Panel (a) shows the logarithmic negativities Em of an initial 
mixed two-mode squeezed state (solid line), a two-mode state 
after local photon subtractions (dotted line) and a two-mode 
state after local photon subtractions combined with optimal 
local displacements (dashed line). Panel (b) shows the success 
rates P SUC c of the concentration protocol with (dashed line) 
and without (dotted line) displacements. Panel (c) shows 
the optimal displacement a opt . The other parameters read 
A = 0.25 and R = 0.1. 



tions. A simple analysis in Fock basis provided an insight 
into the working principle of the squeezing-enhanced en- 
tanglement concentration scheme proposed by Zhang and 
van Loock [231. We have next shown that a similar en- 
tanglement enhancement can be achieved by employing 
local displacements instead of single- mode squeezing. Re- 
markably, the entanglement is enhanced by a destructive 
quantum interference that suppresses or even completely 
eliminates the otherwise dominant vacuum term. The 
displacement-based scheme is much easier to implement 
experimentally than the scheme with local squeezers al- 
beit at the cost of a somewhat reduced success rate of 
the protocol. Note also that while the entanglement 
concentration schemes involving local displacement or 
local squeezing can produce similar amount of entan- 
glement, the internal structure of the entangled states 
is rather different. In the weak-squeezing limit and up 
to local Gaussian operations the displacement enhanced 
scheme essentially yields an entangled single-photon state 
-4=(|10) — 101)) while the scheme with local squeezers pro- 
duces a state exhibiting almost perfect photon number 
correlations -^((OO) + |11)). 

The analytical and numerical results indicate that the 
entanglement enhancement achieved by single photon 
subtractions combined with optimal local displacements 



is of the order of 1 e-bit. Here the entanglement enhance- 
ment is defined as a difference between entanglement of 
concentrated and input states. Even higher entangle- 
ment enhancement could be achieved by more sophisti- 
cated local filtering operations such as the recently pro- 
posed probabilistic noiseless amplifier that approximates 
the unphysical operation g n , where g > 1 is the ampli- 
fication gain [18H24J . Such amplifier may in principle 
generate arbitrarily large amount of entanglement from 
any initial pure two-mode squeezed vacuum state. How- 
ever, the resulting scheme may become rather complex 
and its success rate very small. 

From a more general perspective, the entanglement 
concentration protocol analyzed in this paper is an ex- 
ample of a CV quantum information processing scheme 
whose performance can be significantly enhanced by a 
suitable preprocessing with Gaussian operations. An- 
other example of such improvement is the recently pro- 
posed interface between two weakly coupled quantum 
systems [3l|, |32[ where a perfect transmission of quan- 
tum state from one system to the other can be achieved 
with the help of ancilla modes and auxiliary Gaussian 
operations. It will be interesting to investigate whether 
similar Gaussian pre- or post-processing may be benefi- 
cial also in other CV QIP protocols. 

Note added: After this work was completed we became 
aware of a related paper by Lee et al. [331 ] where entan- 
glement concentration of two-mode squeezed vacuum by 
linear combination of photon addition and subtraction is 
investigated. This protocol is equivalent to the proto- 
col discussed in Sec. II of the present paper, and it is 
shown in Ref. [33] that one ebit of entanglement can be 
obtained by this method in the limit of weak squeezing. 



Acknowledgments 

We acknowledge the financial support of the Future 
and Emerging Technologies (FET) programme within the 
Seventh Framework Programme for Research of the Eu- 
ropean Commission, under the FET-Open grant agree- 
ment COMPAS, number 212008, co-financed by the 
Czech Ministry of Education (7E08028). This work 
was also supported by the Czech Ministry of Edu- 
cation under the projects Center of Modern Optics 
(LC06007) and Measurement and Information in Optics 
(MSM6198959213). 



[1] J. Eisert, S. Scheel, and M.B. Plenio, Phys. Rev. Lett. 

89, 137903 (2002). 
[2] G. Giedke and J.I. Cirac, Phys. Rev. A 66, 032316 

(2002). 
[3] J. Fiurasek, Phys. Rev. Lett. 89, 137904 (2002). 
[4] L.M. Duan, G. Giedke, J.I. Cirac, and P. Zoller, Phys. 

Rev. Lett. 84, 4002 (2000). 
[5] T. Opatrny, G. Kurizki, and D.-G. Welsch, Phys. Rev. 



A 61, 032302 (2000). 
[6] P. T. Cochrane, T. C. Ralph, and G. J. Milburn, Phys. 

Rev. A 65, 062306 (2002). 
[7] S. Olivares, M. G. A. Paris, and R. Bonifacio, Phys. Rev. 

A 67, 032314 (2003). 
[8] A. Ourjoumtsev, A. Dantan, R. Tualle-Brouri, and P. 

Grangier, Phys. Rev. Lett. 98, 030502 (2007). 
[9] H. Takahashi, J.S. Neergaard-Nielsen, M. Takeuchi, M. 



Takeoka, K. Hayasaka, A. Furusawa, and M. Sasaki, Na- [21] 
ture Photonics 4, 178 (2010). 

[10] S.L. Zhang and P. van Loock, Phys. Rev. A 82, 062316 [22] 
(2010). 

[11] A. Ourjoumtsev, R. Tualle-Brouri, J. Laurat, and Ph. 

Grangier, Science 312, 83 (2006). [23] 

[12] J.S. Neergaard-Nielsen, B.M. Nielsen, C. Hettich, K. 
Molmer, and E.S. Polzik, Phys. Rev. Lett. 97, 083604 
(2006). [24] 

[13] K. Wakui, H. Takahashi, A. Furusawa, and M. Sasaki, 

Opt. Express 15, 3568 (2007). [25] 

[14] T. Gerrits, S. Glancy, T.S. Clement, B. Calkins, A.E. [26] 

Lita, A.J. Miller, A.L. Migdall, S.W. Nam, R.P. Mirin, 
and E. Knill, Phys. Rev. A 82, 031802(R) (2010). [27] 

[15] J. Fiurasek, R. Garcfa-Patron, and N.J. Cerf, Phys. Rev. 
A 72, 033822 (2005). 

[16] J.S. Neergaard-Nielsen, M. Takeuchi, K. Wakui, H. Taka- [28] 
hashi, K. Hayasaka, M. Takeoka, and M. Sasaki, Phys. 
Rev. Lett. 105, 053602 (2010). [29] 

[17] P. Marek and J. Fiurasek, Phys. Rev. A 82, 014304 

(2010). [30] 

[18] T.C. Ralph and A.P. Lund, in: Proc. of the Ninth Inter- 
national Conference on Quantum Communication, Mea- [31] 
surement, and Computing, Calgary, Canada, August 9- [32] 
24, 2008 (QCMC 2008), p. 155. [33] 

[19] P. Marek and R. Filip, Phys. Rev. A 81, 022302 (2010). 

[20] J. Fiurasek, Phys. Rev. A 80, 053822 (2009). 



G.Y. Xiang, T.C. Ralph, A.P. Lund, N. Walk, and G.J. 
Pryde, Nature Phot. 4, 316-319 (2010). 

F. Ferreyrol, M. Barbieri, R. Blandino, S. Fossier, R. 
Tualle-Brouri, and P. Grangier, Phys. Rev. Lett. 104, 
123603 (2010). 

M.A. Usuga, C.R. Muller, C. Wittmann, P. Marek, R. 
Filip, C. Marquardt, G. Leuchs, and U.L. Andersen, Na- 
ture Phys. 6, 767 (2010). 
A. Zavatta, J. Fiurasek, and M. Bellini, Nature Phot. 5, 

52 (2011). 

S.L. Zhang and P. van Loock. [arXiv:1103.4500l (2011). 
R. Filip, P. Marek, and U.L. Andersen, Phys. Rev. A 71, 
042308 (2005). 

J.-I. Yoshikawa, T. Hayashi, T. Akiyama, N. Takei, A. 
Huck, U.L. Andersen, and A. Furusawa, Phys. Rev. A 
76, 060301 (R) (2007). 

C.H. Bennett, H.J. Bernstein, S. Popescu, and B. Schu- 
macher, Phys. Rev. A 53, 2046 (1996). 

G. Vidal and R.F. Werner, Phys. Rev. A 65, 032314 
(2002). 

R. Garcfa-Patron, J. Fiurasek, and N.J. Cerf, Phys. Rev. 

A 71, 022105 (2005). 

R. Filip, Phys. Rev. A 80, 022304 (2009). 

P. Marek and R. Filip, Phys. Rev. A 81, 042325 (2010). 

S.-Y. Lee, S.-W. Ji, H.-J. Kim, and H. Nha, Phys. Rev. 

A 84, 012302 (2011). 



